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The rotational invariants constructed by the products ol three spherical harmonic polynomials 
are expressed generally as homogeneous polynomials with respect to the three coordinate vectors, 
where the coefficients are calculated explicitly in this paper. 
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en ' 

r— v I. INTRODUCTION 

^h ' 

Weyl (see p. 53 of [1]) established a theorem on the important structure for rotational invariants: Every even invariant 
~j ■ depending on n vectors ri, r 2 , . . ., r„ in the three-dimensional space IZ3 is expressible in terms of the n 2 scalar products 
Cd ' r a • iv Every odd invariant is a sum of terms 

s ■ 

l_j' [(r xr 6 )-r c ]I(ri,r 2 ,...,r n ), (1) 

t— I [ where r , !■{,, r c are selected from ri, r 2 , ■ ■ ■ , r„, and I(ti, r 2 , ■ • • , r„) is an even invariant. Due to the property of the 
rotational group SO(3), the expression for the invariant J(ri,r 2 , . . . ,r n ) is not unique except for n < 3 because it 
depends upon the coupled orders of n vectors r a . In their famous Encyclopedia of Mathematics on Angular Momentum 
in Quantum Physics [2] , Biedenharn and Louck studied the most important case n = 3 of the general theorem in some 

^s ' detail (§6.17 of [2]), and defined the even invariant (see (6.153) of [2|) as 

O- Wri,r 2 ,r 3 )- / ,4 " 1 ' 



<J2(i k v ! ) ^(riR fc (r 2 )^(r 3 ) (2) 



CM- J " V (2j + l)(2fc+l)(2£+l) 

3 k £ 
s A* v p , 

H ■ (a) a<b 

( j k £\ 

where 1 is the Wigner 3-j symbol, p has to be equal to —(i — v owing to the property of the 3-j symbol, yii?) 

denotes the spherical harmonic polynomial, and ^2 a ab = (j + k + £)j1. The odd invariant Ij,k,£ is proportional to an 
even invariant multiplied by a factor (ri x r 2 ) • r 3 . However, Biedenharn and Louck pointed out in their book (p. 308 of 
2]) that: "Unfortunately, the expression for the general coefficients (6.157) has not been given in the literature, and 
one has had to work out these invariant polynomials from the definition, Eq. (6.153)", where the coefficients (6.157) 
mean Ar a > . 

As a matter of fact, there are only (3rt — 3) independent rotational invariants constructed from n coordinate vectors 
T a . The n(n + l)/2 different scalar products r a ■ r b are obviously too many to be independent for n > 4. The 
independent and complete set of rotational invariants constructed from n coordinate vectors r a is [3| : 

ri -r , r 2 ■ r b , (n x r 2 ) • r c , (3) 
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where 1 < a < n, 2 < b < n, and 3 < c < n. When n = 3, the set © becomes 

?7i=r 2 -r 3 , J72=r 3 -ri, r) 3 = r 1 -r 2 , 

Ci=ri-ri, 6 = r 2 -r 2 , (= (rj xr 2 )-r 3 , (4) 

where C has odd parity and the remaining have even parity. Due to the identity 

C 2 = 666 - 6??i - 6??2 - 6?? 2 + ^mmm, (5) 

with 6 — r 3 ' r 3i the invariant £ may be replaced by 6 in © f° r the even invariant Ij k£{ T ii r 2i r a)' 

The rotational invariant Ij,k,e has widespread applications in many branches of physics, such as in atomic physics [J, 
l5| , in nuclear reaction [2j , in condensed matter physics [6| , in cosmology [7] , and in astronomy and astrophysics |8H10l| . 
Thus, it is desirable to find an explicit expression for Ij.k.e in terms of the rotationally invariant variables rji and 6 

At first sight, the calculation for Ij k,t may be simplified by choosing a special rotation because the spherical 
harmonic polynomial 3^(r) may become much simpler when r points to some special direction, say along the z-axis. 
Making a rotation such that ri points to the z-axis and r 2 is in the xz plane with non- negative x-component, one 
is able to express the invariant Ij,k,£ in terms of the product of special functions (see Appendix). With this idea, 
Harris [11| calculated the even invariants Ij.k.e- His expression for the even Ij.k.e is written in terms of the angles 
among the three vectors ri, r 2 and r3, but not in terms of the scalar products r a • r^ directly. Furthermore, his result 
contains coefficients Cj which are not given explicitly (see (20) in [ll|). Finally, Harris did not calculate the case of 
the odd invariants Ij,k,t- It should be noted that the Harris approach of using spherical variables does not provide the 
coefficients of the invariant polynomials explicitly so that his method does not improve the way of working out these 
invariant polynomials from the definition, as pointed out by Biedenharn and Louck (p. 308 of [2]). 

The purpose of this paper is to present an independent calculation of the coefficients Ai a \ for both even and odd 
invariants Ij.k.e using group theoretical method. The general properties of these invariants are listed in Sec. [TTJ The 
coefficients A( a ) for even and odd invariants are calculated in Sec. IHII and Sec. IIV[ respectively. The conclusions are 
given in Sec. |VJ In Appendix the invariants are expressed in terms of the angles and lengths among three coordinate 
vectors. 

II. GENERAL PROPERTIES OF THE INVARIANTS 

For any given three non-negative integers j, k, and £, satisfying the "triangle rule": 

\j ~ k\ < £ < j + k, (6) 

the rotational invariant Ij.k.ei^ii r 2 ,r3) constructed from the products of three spherical harmonic polynomials is 
defined in ©. The invariant Ij.k.e (ri, r 2 ,rs) has the following properties. 

a) Ij,k,e{Yi, r 2 , v 3 ) is a homogeneous polynomial of orders j, fc, and £ with respect to the coordinate vectors ri, r 2 , 
and r 3 , respectively. 

b) The parity of I h M (n,r 2 ,r 3 ) is (-iy+ k+e . 

c) Due to the symmetry of the Wigner 3-j symbol 



(-1) 



j+k+e 



j k i\ _ ( k j £ 
H v p ) ~ \v n p 



j £ k\ _ ( j k £ 
p p v )~\-ix -v -p 

we have 

(-l)J+fc+^ )M ( ri ,r 2 ,r 3 ) = Ifc^Cra.n.rg) 

= Ij.e. k (ri,r 3 ,r 2 ). (7) 



Thus, we only need to consider Ij,k,e with j < k < i. For the sake of convenience, we write an even invariant as 
Ij,k,j+k~2n and an odd invariant as Ij+i l k+i,j+k-2n+i where < In < j < k. 
d) Ij,k,j+k-2n is real and Ij + i t k+i,j+k-2n+i is pure imaginary because 



E ( H l ) ^(n)^(r 2 )^(r 3 ) 



j k £ 

p v p 
e) Ij ! fe,£(ri,r2,r 3 ) satisfies the three Laplace's equations with respect to ri, r 2 , and r 3 , respectively, 



£ ( J I P ) yu^y-My-oM- 

Hup ^ ' 



A x I jjk/ = A 2 Ij, k / = A 3 7i,M = 0- (8) 

f) From group theory on SO (3) [2J, the decomposition of the direct product of three irreducible representations of 
SO(3), Di{R) x D k (R) x D e (R), where j, fc, and £ satisfy the triangle rule ©, contains one and only one identity 
representation D°(R). Thus, a homogeneous polynomial of orders j, fc, and £ with respect to the coordinate vectors 
ri, r-2, and r 3 , respectively, which satisfies the Laplace's equations (|5J, does exist and is unique up to a constant factor. 

From the explicit formula for the spherical harmonic polynomial Q 



vk , \ k(k + l)(2k+l) ( . ,. . fe _ x , , 
yii(T)=T\/— ^ '-{(x±iy)z h 1 + ...} 



3^(r) = V^^{^ + -}. (») 

where x, y, and z are the three components of r. Among all spherical harmonic polynomials y„(v) with k given, 

the term z k appears only in 3^o( r ) an d the terms xz k ~ 1 and yz k ~ 1 appear only in 3^±i(r)- Therefore, z{z k z^ '~ 

is contained only once in the homogeneous polynomial Ij y k,j+k-2n( r i,^2, T 3), at 3^o( r i)^o ( r 2)3 ; o + _ ™( r 3), with the 
coefficient 



j k j + k — 2n 


Thus the even invariant defined in ([2]) can be rewritten more explicitly in the form 



(10) 



t i \_{jkj + k-2n 

ij,fe,j+fc-2n(.ri,r 2 ,r 3 j = I q q q 

[j/2] j-2a l(k-b)/2] 

X P j,k,j+k-2n 2-^L^, Z^ Aabc 

a=0 6=0 c=max{0,n— a — 6} 

.. /-a/-c/-a+b+c—nk—2c—b') — 2a—bh 



Pj,k,j+k-2n - z2z2z2 A abc, (11) 

a 6 c 

where [m] denotes the largest integer equal to or less than the non-negative real number m. Simi- 
larly, Xiy2z[z2Z^ ~ n is contained only twice in the homogeneous polynomial Ij+i,k+i,j+k-2n+i(^i,^2, r 3 ), at 
yQ 1 (ri)y!g 1 (T2)yt k ~ 2n+1 fa), with the coefficient 



1 (j + 2)\(k + 2)\ ( j + l k + 1 j + k-2n + l , (V)) 



j!fc! V 1 -1 



Thus the odd invariant defined in @ can be rewritten more explicitly in the form 



, v iC /0' + 2)!(fc + 2)l 

Jj+l,fc+l,j+fc-2n+l(,I"l,I"2,l - 3J — -~\ ^T] 

f j + l fc+1 j + k-2n+l 
X { 1 -1 

[j/2] j-2a i(k-b)/2] 

X Qj,k,3+k-2n L^L^i L~i ^ a6c 

a=0 b=0 e=max{0,n-a-h} 

toz-cta+6+c-n fc-2c-6 i-2a—b b 
x Sl?2?3 "l Vl 73; 

Qj,k,j+k-2n = 2^ 2^ 2^ -B a (, c . (13) 

a 6 c 

III. EVEN INVARIANTS 

The coefficients A a b c in (fTTj) are determined from the conditions ((5J, leading directly to the following recursive 
relations: 

2a(2j -2a+ l)A abc 

+ (j -2a-b + 2)(j - 2a - b + l)A (o _ 1)6c 

+(b + 2)(b + l)A {a _ 1){b+2){c _ 1) 

+2(6 + l)(j - 2a - b + l)A {a _ 1)(b+1)c = 0, (14) 

2c(2k-2c+l)A abc 

+ (k-2c-b + 2)(* - 2c - 6 + l)A,6(o-i) 

+(6 + 2)(6 + l)^ (o _i )(w . 2 )(c-i) 

+2(6 + l)(fc - 2c - 6 + l)A»(6+i)(c-i) = 0, (15) 

2 (a + 6 + c - n)(2fc + 2j - 2a - 26 - 2c - 2n + l)A abc 

+(fc - 2c - 6 + 2)(& - 2c - 6 + 1)^60-1) 

+(j - 2a - 6 + 2) (j - 2a - b + l)A {a _ 1)bc 

+2(fc - 2c - 6 + l)(j - 2a - 6 + l)4>(6-i)c = 0. (16) 

Due to the property f) in Sec. UH the solutions for A abc exist uniquely up to a common numerical factor, which can 
be determined for convenience by 

(fc - n)!(2j - 1)!! 



^■OOra — 



(k - 2n)!(2j - 2n - 1)!! 
JT 2m(2j + 2fc-4n-2m + l), (17) 

m— 1 

A = [(j + fc)/2] - n. (18) 



A 

m— 1 



The constant A is chosen to be larger enough so that the coefficients A abc do not contain an unnecessary denominator. 
On the other hand, if A is too large, it may make A abc to be more complex due to the existence of a common factor. 
The choice of A, however, does not effect the final results. 

The coefficients A abc can be calculated one by one by mathematical induction from the recursive relations (I14H16[) . 
where the calculation order is critical. The key point is that at each step of calculation using one of (I14JI16J) . only one 
unknown coefficient is solved from three remaining known coefficients. We choose the calculation order as follows and 
list the calculation results. 

First, we calculate A ob r n _ b \ from (|15[) . and A ab i n _ a _ b \ from (|14[) by mathematical induction: 

{j-n)\{k-n)\ 

*ab(n-a-b) - ^a,b (j _ ^ _ fe) , (fc _ ^ + ^ + &) , 
A 

x ]^[ 2m(2j + 2k - An - 2m + 1), (19) 

m— 1 

where < a < n, < b < n — a, and 

« (_ 1 )q+b+r n! 

' b ~ ^ 2 a - r r\(a - r)\b\ 

(J-2a-b+ r)\(k -2n + 2a + b)\ 

-r)\ 

1MI 
(2j-2n-l)l!(2fc-2n-l)!! ^ 20 ^ 

Evidently, G a ,b = if a < 0, or b < 0, or a + b > n. The function G a ,b will play an essential role for later calculations. 
Second, we calculate A ab i n _ a _ b+C \ from (|16p by mathematical induction: 

b c 

-™-ab{n—a — b+c) / J / J ^a — c+r.b — s 

s— r— max{s,e— a} 

(-l) c 2 s (c!)(j-n)!(fc-n)! 
X (c - r)\{r - s)\s\{j -2a- b)\(k - 2n + 2a + b - 2c)! 

A 

Yl 2m(2j + 2k-An-2m+ 1), (21) 



(n - 2a- b + r)\{j 


-n)\(k-2n + 2a + b-r)\ 


(2j -2a- l)!!(2fc- 


-2n + 4a + 2b-2r-iy.\ 



X 

X 

m— c+1 



where < a < n, < b < n — a, and < c < [(k + b)/2] + a — n. For the case of c > a, terms with c — a > r > s, 
which may occur in the sum over r in (|21[) . vanish because G a ^ b = if a < 0. It is for this reason that the lower bound 
of summation over r in (1211) becomes max{s, c — a}. Similar cases occur in the following formulas. 



Third, we calculate A a r n _ a j rb \ c from (|16p by mathematical induction: 

min{ b-\-c,n~a-{-b} min{ s+c,b+c} 

A a (n-a+b)c = ^ / 4 

s— max{0,6— a} r— max{s,fc+c — a} 

(-l)»+^(b + c)\ 

^ jr a—b—c-\-r.n — a-\-b—s 



(r - s)\s\(b + c - r)\ 

(j — n)\(k — n)\ 
(j — n — a — b)\{k — n + a — b — 2c)! 

A 

x FJ 2m(2j + 2fc-4n-2m + l), (22) 

m— 6+c+l 

where < a < n, < b < j — n — a, and < c < [(k — n + a — b)/2]. For the case of b + c > a, terms with 
s<r<b + c~a, which may occur in the sum over r in (|2"2"j) . vanish because G a ^ — if a < 0. When s > b, terms 
with b + c<r<s + c vanish due to the existence of the factor (b + c — r)! at the denominator. For the case of b > a, 
terms with < s < b — a, which may occur in the sum over s in (|22[) . vanish because G a ,b = if b > n. Finally, for 
the case of c > n — a, terms with n~a + b<s<b + c, which may occur in the sum over s in (|22[) . vanish because 
G a , b = if b < 0. 

Finally, we calculate A( n+0 )j, c from (fT6|) by mathematical induction: 

b s+c 

A( n+a )bc =2^ 2^ 

s— max{0,6— n} r— max{s,b+c— ?i} 

(_l)°+6+c2*( ffl + 6 + c )! 
(r — s)\s\[a + b + c — r)\ 

(j-n)\(k-n)\ 



(j - 2n - 2a - b)\(k - 2c - b)l 

A 

x FJ 2m(2j + 2fc-4n-2m+l), (23) 

m— a+fa+c+l 

where < a < [j/2] — n, < b < j — 2n — 2a, and < c < [(fc — 6)/2]. For the case of 6 + c > n, terms with 
s < r < b + c — n, which may occur in the sum over r in (|23[) , vanish because G Qi ;, = if a < 0. For the case of b > n, 
terms with < s < b — n, which may occur in the sum over s in (|23|) . vanish because G a ,b — if b > n. 



In the following we list some special even invariants Ij.k.j+k-2n for reference 



^0,0,0 — 1, ^0,1,1 — —7= %, 

A),2,2 = ~J= \ ^ [ 3r ?l _ 66] J" , 

^0,3,3 = ^{i[57? 3 -366%]}, 

/o,4,4 = ^ 1 1 [35^ _ 3066r/ 2 + g^2] 

Jo,5,8 = -i (^ [63^ - 7066r?! + l^|»?i] 
Vll 1° 

Jo,6,6 - "i (4 [ 23l7 ?l - 315 66?7l 



03 U6 

+ 105^1^-5^1]}, 



In 



1 123 — 



h 34 — 



'145 



1 

15 12 



[3771772 - 6%] f , 



V 35 1 2 [ 57? ^ 2 " ^ 2 ^ 3?72 ~~ 2 6%%] 

JV 7 I 8 t 35 ^ 772 ~ 15 66%% ~ i5 6?7 2 % 



366 2 %]}, 

1 r5" n 



;;V |;| [6377*772 - 4266%*% 

+ 3£|^2 - 286r? 3 77 3 + 126£f %%] } , 

/l56 = V li3 { 1^ [ 231 ^ 7 ' 2 " 210 66^% + 35£&|%% 
- 1056^773 + 7066V% - 5£&f %] } . 



^222 — — 



2J1 

35 \ 2 
■36% 2 + 2^66] } , 



36»?2 + 9771772773 - 36?7 2 



-^224 = V 35 1 8 t 3577 ^ 2 ~ 5 &&V2 - 206%%% 

+26 2 % 2 -566?7 2 +666 2 ]}, 



^233 — 



ao5 



[-306t7i772 + 75t7 2 ?72773 - 366%% 



-306%% 2 - 256r7 3 + 21666%] } , 
/235 = " V 2^1 \S t 63 "^ " 216e377l7?2 ~ 42 6»7i»&»/3 

+666 2 %% + 66 2 %% 2 - 766??? + 3666 2 %] } . 



' 2 .1. 1 



1 

77 18 



636^1% + 966% + 14777^2773 



- 2766%%% - 636ry 2 77 3 2 + 96£f r?f " ^Zivi 

+ 5i666r??-666 2 6 2 ]}, 



1686%%% + 566^|%%% + 28^77^77^ 

4f„f3 T1 2 _01^f„„4 , ^At.tt2 n 2 _f,t2rfll 



[2 ( 1 

/ 334 = -y 77 { s \r 1Q ^™ + 1757 ? 2? fe 2 % + 566^% 
+60666%% + 566% 2 % - 7666 2 %] } - 

I 3 36 - -^ {^ [231^^3 - 63&&W - 1896^ 2 % 2 % 
+ 2166% 2 % + 42$ViV2Vl - 2^4 - 6366^1% 

+2i666 2 '7i% + 2i66 2 ^ 2 % - 3666 3 %] } • 



IV. ODD INVARIANTS 

Substituting (|T^|) into the Laplace's equations ©, we obtain the recursive relations for the coefficients B abc 

2a{2j -2a + i)B abc 

+ (j -2a-b + 2)0" - 2a - b + l)B (o _ 1)6c 

+(6 + 2)(6+l)B (o _i )(6+2)(o _i ) 

+2(6 + l)(j - 2a - 6 + l)B (a _i )(6+1)c = 0, (24) 

2c(2fc-2c + 3)B abc 

+ (fc - 2c - b + 2)(k - 2c - b + l)-B o6(c _i) 

+(6 + 2)(6 + l)B (o _i )(6+ 2) (e _i) 

+2(6 + l)(fc - 2c - b + l)B (5+i)(c-i) = 0, (25) 

2(a + b + c - n)(2k + 2j - 2a - 2b - 2c - 2n + 3)B abc 
+(k-2c-b + 2){k -2c-b+ l)B ab ( c _i) 
+(j -2a-b + 2)(j-2a-b+ l)B (a _ 1)bc 

+2(fc - 2c - b + 1)(>- - 2a - b + l)5 o(5 _i )c = 0. (26) 

The only difference between the two sets of equations (J14M16I) and (|24j|26[) is that "(+1)" in the first term of each 
equation of the first set is replaced by "(+3)" in that of the second set. Thus, through the same procedure we have 
calculated the coefficients B abc listed below, which are very similar to the coefficients A abc . 

_ (fc-n)!(2j + l)!! 

-Dnnri 



*3ab(n—a — b) * 'a, 



(fc-2n)!(2j-2n + l)!! 

J[ 2m(2j + 2fc - An - 2m + 3). (27) 

(j — n)\(k — n)\ 



x 

X 

771—1 



(j-2a-b)\{k-2n + 2a + b)\ 
Y[ 2m(2j + 2k-in-2m + 3), (28) 



A 
X 

777 — 1 



where < a < n, < 6 < n — a, and 

(-l) a+b+r n! 



Fa, b = £ 



2 a - r r!(a-r)!6! 

r=max{0 : 2a+6-n} v ' 

(j -2a-b + r)\(k -2n + 2a + b)\ 
(n-2a-b + r)l(j - n)\(k -2n + 2a + b-r)\ 

(2j -2a + l)!!(2fc - 2n + 4a + 2b- 2r + 1)!! 

X (2j - 2n + l)\\(2k - 2n + l)\\ ' 

Evidently, F a ^ b = if a < 0, or b < 0, or a + b > n. 

b c 

£>ab(n — a — b+c) = / d / J & a — c+r,b—s 

s— r— max{s,c— a} 

(-l) c 2 s (c!)(j-n)!(£;-n)! 



A 

X 

m— c+1 



where < a < n, 0<b<n-a, and < c < [(fc + b)/2] + a 



n. 



min{6-fc,n — a+6} min{s+c,6+c} 

S (n-«+6)c = 2^ / y 

s— max{0,6— a} r— max{s,6+c— a} 

(-l) h+c 2 s (6 + c)! 

"a — b— c-\-r.n — a-\-b— s 



where 0<a<n, 0<6<j — n — a, and < c < [(& — n + a — b)/2]. 

b s-\-c 

B( n +a)bc 22 22 

s— max{0,6- n} r-max{s,Hc-n} 

(-l) Q+fc+c 2 s (a + 6 + c)! 

F n — b— c+r,b— s 



(r — s)\s\(a + b + c — r)! 
C? — w)!(fc — n)! 



A 

x 



(29) 



(c - r)!(r - s)\s\(j -2a- b)\(k -2n + 2a + b- 2c)! 

A 

rj 2m(2j + 2k-4n-2m + 3), (30) 



(r - s)!s!(6 + c-r)! 

(.7 — w)!(fe — n)\ 

(j — n — a — b)\(k — n + a — b — 2c)! 

A 
x J| 2m(2j + 2fc-4n-2m + 3), (31) 

?n— 6+c+l 



(j -2n- 2a -&)!(& -2c -6)! 

A 

J| 2m(2j + 2fc-4n-2m + 3), (32) 



10 



where < a < [j/2] - n, < b < j - In - 2a, and < c < [(k - b)/2\. 

In the following we list some special odd invariants Ij+\,k+\ ,j+k-2n+i for reference. 



11 



^1,1,1 

Jl,2,2 
h, 3, 3 

^1,4,4 
h, 5, 5 
-^1,6,6 

h,7,7 

^2,2,3 
^2,3,4 
-^2,4,5 

-^2,5,6 

-^2,6,7 

-^3,3,3 

^3,3,5 



3] r • 



= - iC Vio ^' 

= icy^fl || [33^ - 3066% 3 + Hkhi] 

= iC \/l{6i [429?? ^ 49566r/ ' 
+ 135C 2 2 C 3 V-5^3 3 ]}, 

= iC3W — <^ - [SviVi - 6%] J- , 

= -iCy -j \ j [?vlv2 - 66% - 2^3771773] 

/TO firs 2 

= iCy YT J 8 L 21rM2 ~ 7 66%% ~ 7 &%% 
+ 66 2 %]}, 

= -&J§£ [l t 33 ^ - 1866^% 

+ elel/72 - I26r??% + 466 2 %%] } , 



= 2-^11 {-^ [429^2 - 33066^?% + 456^%% 
- 165^^% + 9066V% - 5^1 %] } , 

= -iCy j \ y^ [- 5 6»?2 + 25r/i 773773 - 56?? 2 
-56 77? + 2666] } , 

= iC5^ j-^ [63r? 2 r? 2 - 766% 2 - 286771772773 
+26 2 % 2 -766r/ 2 +666 2 ]}, 

- iC3^ |^ [-706771 t? 2 + 245772772% 

-1566%% - 706t?i 77! - 496t?? + 31666%] } , 



73,4,6 



'3,5,5 



+266 2 %% + 2£f %r?f - 366% 3 + 666 2 %] } , 

= -^yH{^ [- 45 6r? 2 % 2 + 56 2 6% 2 + 135r7?77 2 r?3 
- 2566%%% - 456t? 2 % 2 + 566 2 % 2 - 276r?f 
+ 23666^7? - 26 ZUZ]}, 
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- 336^^ + mi^hi - Ci^el] } > 

/ 4 , 4 ,5 = -iC^= { 4 [-1266m??2 + 44177^1,73 - 7&&^3 
V143 L'2 

- 1966m»72?7| + 14e|??| - 126Ci7??772 

14 / 10 f 1 r T o , on 

4,4,7 = iCy y Y43 |g4 [4291??^ - 99£ 2 £ 3 W - 2976^% 
+ 276C 3 2 % 2 % + 54$77i 772771 - 2$ »jf - 99&&»tf»fc 

+wtiitehim + mie 3 vh3 - 366^1%] } ■ 



V. CONCLUSIONS 

The rotational invariants Ij,k,t(*ii i"2 ; r 3 ) constructed by three spherical harmonic polynomials are the homogeneous 
polynomials of orders j, k, and £ with respect to the three coordinate vectors r\, r 2 , and r3, respectively. We have 
rewritten the definitions for the invariants given by Biedenharn and Louck more explicitly, derived the recursive 
relations for the coefficients A a b c and B a b c using the Laplace's equations, defined two key functions G a ,b and F a ^ in 
(|20|) and (|29j) . and calculated analytically the expressions for A a \, c and B a i, c by mathematical induction, as given in 
(J2HI23J) and (J30M32J) . Therefore, we have completely solved the problem raised by Biedenharn and Louck (p. 308 of [2]). 
The present method can in principle be generalized to the rotational invariants constructed by four or more spherical 
harmonic polynomials although the definition for the invariants depends on the order of coupling. 

Acknowledgments 

ZQM would like to thank Professor Fu-Chun Zhang for his warm hospitality during his visit at the University of 
Hong Kong where part of this work was completed. This work was partly supported by the Natural Science Foundation 
of China under grant No. 11174099 and No.11075014. ZCY was supported by NSERC of Canada. 



[1] H. Weyl, The Classical Groups: Their Invariants and Representations, (Princeton Univ. Press, Princeton, N. J. 1946). 

[2] L. C. Biedenharn and J. D. Louck, Angular Momentum in Quantum Physics, Theory and Application, Encyclopedia of 
Mathematics and its Application, Ed. G.-C. Rota, (Addison- Wesley. Massachusetts, 1981). 

[3] X.-Y. Gu, B. Duan and Z.-Q. Ma, Phys. Rev. A 64, 042108 (2001). 

[4] D. M. Fromm and R. N. Hill, Phys. Rev. A 36, 1013 (1987). 

[5] N. L. Manakov, A. V. Meremianin, and A. F. Starace, Phys. Rev. A 57, 3233 (1998). 

[6] O. Borisenko and V. Kushnir, Ukr. J. Phys. 51, N 1, 90 (2006). 

[7] N. Joshi, S. Jhingan, T. Souradeep, and A. Hajian, Phys. Rev. D 81, 083012 (2010). 

[8] G. Rocha, M. P. Hobson, S. Smith, P. Ferreira, and A. Challinor, Mon. Not. R. Astron. Soc. 357, 1 (2005). 

[9] P. Papai and I. Szapudi, Mon. Not. R. Astron. Soc. 389, 292 (2008). 
[10] A. Raccanelli, L. Samushia, and W. J. Percival, Mon. Not. R. Astron. Soc. 409, 1525 (2010). 
[11] F. E. Harris, Integrals for exponentially correlated four-body systems of general angular symmetry, in Fundamental World 

of Quantum Chemistry, Vol. 3, p. 115, Ed. E. J. Brandas and E. S. Kryachko, (Kluwer, Dordrecht, 2004). 
[12] I. S. Gradshteyn, I. M. Ryzhik, Ta67e of Integrals, Series, and Products, 7th Ed., Editors: A. Jeffrey and D. Zwillinger 
(Academic Press, 2007). 



13 



Appendix: Invariants expressed by spherical variables 

Choosing a special rotation such that ri is along the z axis and v 2 is in the xz plane with x positive, we have 

?(rO=^(d /2 ,0 ) 0) = /^I^ ) 



y 

yi(v 2 ) = yt(i l 2'\o l2 ,Q) 

= r _ n( v +M )/2 [ (2A: + l)e 2 fc (fc-H)!(fc + H)! 
1 ; 47r4l"l(fc!) 2 



x (sinfli 2 )l"lpWjW(cosfli 2 ) 



1/2 



= r_n(-+i-i)/a r (2^ + i)^-H)!(^ + H)! l] J 

4tt4IH(^!)2 



x(sin0 13 )^lpl^ l (cos0 13 )e 
where Pn is the Jacobi's polynomial (see 8.960 in [12 ]) 



-11/^3 



and the angles satisfy 



Thus, 



m=0 



ii + a 
m 



m=0 x ' x 

(x -!)"-"> (x + 1) 



n + /3 
n — m 



3 

COS ab = ^ ZabcVc {iaibT 1 , 
c=\ 
COS #23 — COS 012 COS 013 

cost/? = — — , 

sin 0i2 sin 0i3 

C = (Ci^26) 1/2 sin6» 12 sin6'i3sin93. 



From the identity 



Ij,kj{ri,r 2 ,r 3 ) = 



zkki ^ 



kW. 






V(fc - i/)!(fc + v)\(e - v)\(i + ljy. , 

4l"l ° 



sin 0i2 sin 13 ) ' ' 



l-\ v \ (cos 6 12 )Pj,^ (cos 9 



i3je 



— ii/tp 



COS 



(vip) + i sin (vip) = e ll " p = [cos ip + i sin <p] v 
= H ( r ) ( cos v) !/ " r ( isin ^) r , 



(A.1) 



14 
one obtains for v > 

cos*, -eb-W^.) (^)(cos^r 2r+2s , 

r=0 s=0 ^ / \ / 

[("-l)/2] r , s 

sini/p = sin^ ^ I^" 1 )^ ( 2r + 1 ) 



r \ (n nB ,A"- 1 -» r + 2 ' 



Substituting them into (A.l), wc obtain 



X | g I (COS^) 



for even j + fc + t, and 



for odd j + k 



J, iM (r 1 ,r 2 ,r 3 ) = -^ r - £ (-1) (^ J ^ _„ J 

x jh p k-\l\(cose 12 ) 

[1-1/2] r / I I \ / \ (A ' 2) 



r=0 s=0 
X (COS #23 — COS012 COS 6*13) 



[(l-cos 2 fl 12 )(l-cos 2 13 )p s , 

^■,fe^(ri,i-2,r 3 ) = -i£-¥— i — 1| — ^ (-l)"sign(z/) 

/ 3 k i \ V \k - i/)!(fc + gig - i/)!(l +g! 
y zv -1/ y 4IH 



i_\l\ (cos6>i 2 )P;_^| l (cos(9i3) 



r fc-M ^ Uh(7 12^£-|„| ^w»"i3; (A.3) 

I(k|-l)/2] r , 

r=0 s=0 v ' v 

X (COS 6*23 — COS 0\2 COS 0\3) 



[(l-CO S 2 ^2)(l-COS 2 ^3)] r " S . 



